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Most of the analytic approaches which are used at present to understand the low energy hadronic interactions 
in Particle Physics, get their inspiration from QCD in the limit of a large number of colors N c . I first illustrate 
this with the example of the left-right correlation function which is an excellent theoretical laboratory. Next, I 
present the list of observables which have been computed using a large-iV c QCD approach. Finally, I discuss in 
some detail examples which are relevant to lattice QCD, in the sense that we can make comparisons. 



1. INTRODUCTION 

In the Standard Model, the electroweak in- 
teractions of hadrons at very low energies are 
conveniently described by an effective chiral La- 
grangian, which has as active degrees of freedom 
the low lying SU(3) octet of pseudoscalar parti- 
cles, plus leptons and photons. The underlying 
theory is a SU(3) C x SU(2) L x U(l) w gauge the- 
ory which is formulated in terms of quarks, glu- 
ons and leptons, together with the massive gauge 
fields of the electroweak interactions and the hith- 
erto unobserved Higgs particle. Going from the 
underlying Lagrangian to the effective chiral La- 
grangian is a typical renormalization group prob- 
lem. It has been possible to integrate the heavy 
degrees of freedom of the underlying theory, in 
the presence of the strong interactions, pertur- 
batively, thanks to the asymptotic freedom prop- 
erty of the £t/(3)-QCD sector of the theory. This 
brings us down to an effective field theory which 
consists of the QCD Lagrangian with the u, d, s 
quarks still active, plus a string of four quark op- 
erators and mixed quark-lepton operators, mod- 
ulated by coefficients which are functions of the 
masses of the fields which have been integrated 
out and the scale /x of whatever renormalization 
scheme has been used to carry out this integra- 
tion. We are still left with the evolution from 
this effective field theory, appropriate at inter- 
mediate scales of the order of a few GeV, down 
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to an effective Lagrangian description in terms 
of the low-lying pseudoscalar particles which are 
the Goldstone modes associated with the spon- 
taneous symmetry breaking of chiral-5'J7(3) in 
the light quark sector. The dynamical descrip- 
tion of this evolution involves genuinely non- 
perturbative phenomena and it is mostly studied 
using the techniques of Lattice QCD. In this talk, 
I shall review the progress which has been made 
in approaching this last step using analytic meth- 
ods; in particular when the problem is formulated 
within the context of QCD in the limit of a large 
number of colours N c . 

The suggestion to keep N c in QCD as a free 
parameter was first made by G. 't Hooft ]jj as 
a possible way to approach the study of non- 
perturbative aspects of QCD. The limit N c — > 
oo is taken with the product a s N c kept fixed. 
In spite of the efforts of many illustrious the- 
orists who have worked on the subject, QCD 
in the large- N c limit (QCD^) still remains un- 
solved; but many interesting properties have 
been proved, which suggest that, indeed, the 
theory in this limit has the bulk of the non- 
perturbative behaviour of full QCD. In particular, 
it has been shown that, if confinement persists in 
this limit, there is spontaneous chiral symmetry 
breaking 0] . It should be stressed that QCD^ is 
not just a "wild extrapolation" from N c = 3 to 
N c — oo. In fact, N c is really used as a label 
to select specific topologies among Feynman di- 
agrams. The topology which corresponds to the 
highest power in the iV c -label is the one which se- 



2 



lects planar diagrams only; and the claim is that 
this class provides already a good approximation 
to the full theory. 

The spectrum of QCD^ consists of an infinite 
number of narrow stable meson states which are 
flavour nonets j|. This spectrum looks a pri- 
ori rather different to the one in the real world: 
the examples of the vector and axial-vector spec- 
tral functions measured in e + e~ — > hadrons and 
in the hadronic r-decay which are shown below, 
have indeed a richer structure than just a sum 
of narrow states. There are, however, many in- 
stances where one can show that the observables 
one is interested in, are given by smooth integrals 
of specific hadronic spectral functions, though un- 
fortunately in most cases, the specific hadronic 
spectral functions in question are not accessible 
from experiments. Typical examples of such ob- 
servables, are the coupling constants of the effec- 
tive chiral Lagrangian of QCD at low energies, 
as well as the coupling constants of the effective 
chiral Lagrangian of the electroweak interactions 
of pseudoscalar particles in the Standard Model. 
What is needed for the evaluation of these cou- 
pling constants is not so much the detailed point- 
by-point knowledge of the relevant hadronic spec- 
trum, but rather a good approximation consistent 
with the asymptotic properties of QCD, both at 
short- and long-distances. It is in this sense that 
the simple QCDgo-spectrum becomes useful. It 
provides a simple parameterization of the physi- 
cal hadronic spectrum, based on first principles. 

2. THE CHIRAL LAGRANGIAN 

The strong and electroweak interactions of the 
Goldstone modes at very low energies are de- 
scribed by an effective Lagrangian which has 
terms with an increasing number of derivatives 
(and quark masses if explicit chiral symmetry 
breaking is taken into account.) Typical terms 
of the chiral Lagrangian are 

C eS = ~F 2 tr (DvUIVUt) (1) 

V v ' 

7T7T — >mr , K — tireis 

+ L 10 tr(U^F Rflv UF^)+--- (2) 



+ e 2 Ctr(Q R UQ L U^) +■■• (3) 

S ' 

- e 2 C -^(TT+TT-+K+K-) 

- ^V ud V: s g s _F* (D^UD^)^ + ■■■, (4) 

K — >7T7T , K — >7T7r7r 

where U is a 3 x 3 unitary matrix in flavour space 
which collects the Goldstone fields and which un- 
der chiral rotations transforms as U — > V R UV£; 
D^U denotes the covariant derivative in the pres- 
ence of external vector and axial-vector sources. 
The first line is the lowest order term in the 
sector of the strong interactions [Q, Fo is the 
pion-decay coupling constant in the chiral limit 
where the light quark masses u, d, s are neglected 
(Fo ~ 90 MeV); the second line shows one of 
the couplings at 0(p 4 ) |^,^); the third line shows 
the lowest order term which appears when pho- 
tons and Z's are integrated out (Ql = Qr = 
diag.[(2/3, — 1/3, — 1/3]), in the presence of the 
strong interactions ; the fourth line shows one of 
the lowest order terms in the sector of the weak 
interactions. The typical physical processes to 
which each term contributes are indicated under 
the braces. Each term is modulated by a coupling 
constant: F5 2 , Liq,... C...g%..., which encodes the 
underlying dynamics responsible for the appear- 
ance of the corresponding effective term. The 
evaluation of these couplings from the underlying 
theory is the question we are interested in. The 
coupling #8 for example, governs the strength of 
the dominant AI — 1/2 transitions for F-decays 
to leading order in the chiral expansion. 

2.1. Two crucial observations 

There are two crucial observations to be made 
concerning the relation of these low energy con- 
stants to the underlying theory. 

i) The low-energy constants of the Strong La- 
grangian, like Fg and L\q, are the coefficients 
of the Taylor expansion of appropriate QCD 
Green's Functions. For example, with Hlr.(Q 2 ) 
the correlation function of a left-current with a 
right-current in the chiral limit: 

J d A x e^ x {Q\T{u L ^d L {x)u Rl v d R {0^) |0) 
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= ^(gV-<T<7 2 )n M (Q 2 ), (5) 

the Taylor expansion 

- Q 2 IWQ 2 )| Q ^ = Fl - AL W Q 2 + ■ ■ ■ , (6) 

defines the constants Fq and L±o. 

ii) By contrast, the low-energy constants of the 
Electroweak Lagrangian, like e.g. C and g$, are 
integrals of appropriate QCD Green's Functions. 
For example Q], 

r^Aj- Q ' u ^ Qi)) - <7) 

Their evaluation appears to be, a priori, quite a 
formidable task because they require the knowl- 
edge of Green's functions at all values of the 
euclidean momenta; i.e. they require a precise 
matching of the short-distance and the long- 
distance contributions to the underlying Green's 
functions. 

These two observations are generic in the case 
of the Standard Model, independently of the 
l/iV c -expansion. The large-iV c approximation 
helps, however, because it restricts the analytic 
structure of the Green's functions in general, and 
Hlr(Q 2 ) in particular, to be meromorphic func- 
tions: they only have poles as singularities; e.g., 
in QCD^, 

TT (rP\—X^ fv M V \ " fl M A F Q fax 

ULRiQ 

where the sums are extended to an infinite num- 
ber of states. In practice, however, in the case of 
Green's functions which are order parameters of 
spontaneous chiral symmetry breaking (like U LR 
in the chiral limit), these sums will be restricted 
to a finite number of states. 

2.2. Sum Rules 

There are two types of important restrictions 
on Green's functions like Ulr{Q 2 )- One type 
follows from the fact that, as already stated 
above, the Taylor expansion at low euclidean mo- 
menta must match the low energy constants of 
the strong chiral Lagrangian. This results in a 
series of long-distance sum rules like e.g. 

£#-£# = -^10. o) 

V A 



Another type of constraints follows from 
the short-distance properties of the underlying 
Green's functions. The behaviour at large eu- 
clidean momenta of the Green's functions which 
govern the low energy constants of the chiral La- 
grangian can be obtained from the operator prod- 
uct expansion (OPE) of local currents in QCD. In 
the large- N c limit, this results in a series of al- 
gebraic sum rules || which restrict the coupling- 
constants and masses of the hadronic poles. In 
the case of the Li?-correlation function in Eq. (||) 
one has e.g., 

E/^-E/lMi-F 2 =0, (10) 
J2tiM*-J2flMi = 0, (11) 

£ ftM$ - £ f A M% ~ -4™, (^V) 2 ■ (12) 

The sum rules in Eqs. (|To| ) and ([ll]) are the cele- 
brated Weinberg sum rules, which follow from the 
fact that in the chiral limit, there are no 0(l/Q 2 ) 
terms and no 0(1/Q 4 ) terms in the OPE of 
IWQ 2 ). The third sum rule in Eq. @ §, 
follows from the matching between the 0(1/(2°)- 
terms in the QCD^ expression in Eq. (||) and the 
corresponding one in the OPE (evaluated at lead- 
ing 0(a s N 2 ) M). In principle there are an infi- 
nite number of sum rules in QCD^, which relate 
the masses and couplings of the QCD OT narrow 
states to the local order parameters which appear 
in the OPE and to the non-local order parameters 
which govern the chiral expansion in the strong 
sector. 

3. THE 7T+ tt° MASS DIFFERENCE AS 
A THEORETICAL LABORATORY 

The physical effect of the coupling C in Eq. (|^) 
is to give a mass, mostly of electromagnetic origin, 
to the charged pions: 

™Ulu = ^ d<2 2 [-Q 2 TIlr(Q 2 )] ■ (13) 

The first evaluation of this integral was done by 
F. Low and collaborators in 1967 @. What I 
am going to do next is, simply, to put their phe- 
nomenological calculation within the context of 
QCD ro . 
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As already stated, the function Ulr{Q 2 ) in 
QCD^ is a meromorphic function. That means 
that within a finite radius in the complex Q 2 - 
plane it only has a finite number of poles. The 
natural question which arises is: what is the min- 
imal number of poles required to satisfy the OPE 
constraints? The answer to that follows from a 
well known theorem in analysis jll| , which when 
applied to our case, states that the function 

- Q 2 U LR (Q 2 ) ee A[z] with z = ^L, (14) 

My 

and My the mass of the lowest narrow state, has 
the property that 

J_ f A'\z] 
2ttI 



Af-V 



-dz . 



(15) 



where Af is the number of zeros and V is the num- 
ber of poles inside the integration contour in the 
complex z-plane (a zero and/or a pole of order m 
is counted m times.) For a circle of radius suffi- 
ciently large, we simply have that Af — V = —p , 
where p denotes the asymptotic power fall-off in 
1/z of the A[z] function. Identifying the power 
p with the leading power predicted by the OPE, 
and from the fact that Af > 0, there follows that 
V > Pope- In our case Pope = 2 =>• V > 2 and the 
minimal hadronic approximation (MHA) compat- 
ible with the OPE requires two poles 0: one vector 
state and one axial-vector state. The MHA to the 
QCD^ expression in Eq. (@) is then the simple 
function 



Q 2 IWQ 2 ) - F< 



M 2 M\ 



(Q 2 + M V )(Q 2 + MX) 



.(16) 



Inserting this function in Eq. (O) gives a predic- 



tion for the electromagnetic ir + 
difference, with the result ^ 

Am T = (4.9±0.4)MeV, 



-7T ee Am T mass 



(17) 



to be compared with the experimental value pj ] 
Am^ = (4.5936 ± 0.0005) MeV . (18) 



2 Notice that with the definition of A[z] in Eq. (|14|) the 
pion pole is removed. 

(87 ± 3.5) MeV, My = 
0.50 ± 0.06. These values 
of the low energy 



3 This is the result for Fo 
(748 ± 29) MeV and g A = ^ 



follow from an overall fit to predicti 
constants in the MHA to QCD^, fll 



Figure 1. The Spectral Function -l-ImIT£/j(i), ob- 
tained from the ALEPH data |7j|/ compared to the 
MHA to QCD^. 
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3.1. The real world and the MHA to QCD^ 
The spectral function of the correlation func- 
tion defined in Eq. (|^) can be obtained from 
measurements of the hadronic r-decay spectrum 
(vector-like decays minus axial-vector like de- 
cays). We plot in Figure [l] the experimental 
determination of ilmIT£fl(t), obtained from the 
ALEPH collaboration data (l5) , versus the invari- 
ant hadronic mass squared t in the accessible re- 
gion < t < m 2 . In this case, the correspond- 
ing spectrum of the MHA to QCD^ consists of 
the pion pole (not shown in the figure), a vector 
narrow state(the p) and an axial-vector narrow 
state (the A\). At this level of approximation, 
and in the chiral limit, the rest of the vector and 
axial-vector states are taken to be degenerate and 
cancel in the difference, a fact which is simulated 
by the horizontal continuum line in the figure. 
Looking at this plot, one can hardly claim that 
this approximation reproduces the details of the 
experimental data. However , with Hlr(Q 2 ) de- 
termined from the spectral function by the un- 
subtracted dispersion relation 



IWQ 2 ) 



dt- 



1 1 



, t + Q 2 ir 

the corresponding plot 
versus 



~TmH LR (t) , 



(19) 



of the function 
~W^^lr(Q 2 ) versus the euclidean variable 
Q 2 = — i is shown in Fig g. The solid curve 



5 



is the one corresponding to the simple MHA in 
Eq. d||) and the dotted curve the one from the 
experimental data in Fig [|. One can see that, 
by contrast to what happens in the Minkowski 
region shown in Fig. [I], the corresponding curves 
in the euclidean region are both very smooth and 
in fact the MHA already provides a rather good 
interpolation between the asymptotic regimes 
where, by construction, it has been constrained 
to satisfy the lowest order chiral behaviour in 
Eq. and the two Weinberg sum rules: the 
OPE constraints in Eqs. © and This 
good interpolation is the reason why the inte- 
gral in Eq. (0), evaluated at the MHA, already 
reproduces the experimental result rather well . 

3.2. Other Analytic Approaches 

At this stage, it is illustrative to compare the 
large-iVc approach we have discussed so far with 
other analytic approaches in the literature. The 
TLlr correlation function provides us with an ex- 
cellent theoretical laboratory to do this compar- 
ison. The different shapes of z: M-IiLR.(Q 2 ) pre- 
dieted by other analytic approaches are collected 
in Fig. [| Here, several comments are in order. 

a) The suggestion to use a large--/V c QCD 
framework combined with xPT cut-off loops, was 
first proposed by Bardeen, Buras and Gerard 
in a series of seminal papers [[l6]jl7],[l8| . The 
same method has been applied by the Dort- 
mund group 1 19 1, in particular to the evaluation 
of e'/e. In this approach the hadronic ansatz 
to the Green's functions consists of Goldstone 
poles alone and, therefore, integrals like the one 
in Eq. (^) become UV-divergent (often quadrat- 
ically divergent) since the correct QCD short- 
distance behaviour is not implemented. In prac- 
tice the integrals are cut at a physical cut off: 
A ~ 1 GeV. In the case we are considering, the 
predicted shape of the LR-correlation function 
normalized to its value at Q 2 = is a constant all 
the way up to the chosen cut-off value, as shown 
by the BBG, HKPSB line (the horizontal dotted 
line) in Fig. [|. 

b) The Trieste group evaluate the relevant 
Green's functions using the constituent chiral 
quark model (CxQM) proposed in refs. |^0) and 
plpa] . They have obtained a long list of predic- 



Figure 2. Plot of —^-^Ilr(Q 2 ) in the euclidean 
region. The solid curve is the one correspond- 
ing to the simple MHA in Eq. fildj ) and the dot- 
ted curve the one from the experimental data in 
Fig 0. The other curves are the predictions of 
various models discussed in the text. 



0.6 



0.4 



0.2 



\\ BBG, HKPSB 

\\ 

I \ 




\ \ 




ALEPH '"•^Nw^ 


ENJL 



tions ^3|, in particular e'/e. The model gives an 
educated first guess of the low-Q 2 behaviour of 
the Green's functions, as one can judge from the 
CxQM-curve (dashed curve) in Fig. |2[ but it fails 
to reproduce the short-distance QCD-behaviour. 
Another objection to this approach is that the 
"natural matching scale" to the short-distance 
behaviour in this model should be ~ 4Mq (Mq 
the constituent quark mass), which is too low to 
be trusted. 

c) The extended Nambu-Jona-Lasinio (ENJL) 
model was developed as an improvement on the 
CxQM, since in a certain way it incorporates the 
vector-like fluctuations of the underlying QCD 
theory which are known to be phenomenologically 
important (see e.g. refs. 0,^5) where other refer- 
ences can also be found) . The model is rather suc- 
cessful in predicting the low-energy 0(p 4 ) con- 
stants of the chiral Lagrangian [24f| . It has, in- 
deed, a better low-energy behaviour than the 
CxQM, as the ENJL-curve (dot -dashed) in Fig. | 
shows; but it fails to reproduce the short-distance 
behaviour of the OPE in QCD. Arguments, how- 
ever, to do the matching to short-distance QCD 
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have been forcefully elaborated in refs. [|26|, which 
also give a lot of numerical predictions; a large 
value for e'/e in particular. 

The problem with the ENJL model as a plausi- 
ble model of large-iV c QCD, is that the on-shcll 
production of unconfined constituent quark QQ 
pairs that it predicts, violates the large- N c QCD 
counting rules. In fact, as shown in ref. when 
the unconfining pieces in the ENJL spectral func- 
tions are removed by adding an appropriate series 
of local counterterms, the resulting theory is en- 
tirely equivalent to an effective chiral meson the- 
ory with three narrow states V, A and S; very 
similar to the phenomenological Resonance Chi- 
ral Lagrangians proposed in refs. p7|]2§[ ]. These 
Lagrangians can be viewed, therefore, as partic- 
ular models of large- N c QCD. They predict the 
same Green's functions as the MHA to QCD^ dis- 
cussed above, in some particular cases but not in 
general (see e .g. the three-point functions dis- 
cussed in ref. |29(| ). 

In view of the difficulties which the analytic 
approaches discussed in a), b) and c) above, have 
in reproducing the shape of the simplest Green's 
function one can think of, it is difficult to at- 
tribute more than a qualitative significance to 
their "predictions" ; e'/e in particular, which re- 
quires the interplay of several other Green's func- 
tions much more complex than Hlb.(Q 2 )- 

4. METHODOLOGY, APPLICATIONS 

The QCDqq approach that we are proposing in 
order to compute a specific coupling of the chiral 
electroweak Lagrangian consists of the following 
steps: 

1. Identify the relevant Green's functions. 

In most cases of interest, the underlying 
QCD Green's functions in question are two- 
point functions with additional zero mo- 
mentum insertions of vector, axial vec- 
tor, scalar and pseudoscalar currents. The 
higher the power in the chiral expansion, 
the higher will be the number of zero mo- 
mentum insertions. This step is totally gen- 
eral and does not invoke any large-iV c ap- 
proximation. 



2. Work out the short-distance behaviour and 
the long-distance behaviour of the relevant 
Green's functions. 

The long-distance behaviour is governed by 
the Goldstone singularities and can be ob- 
tained from xPT. The short-distance be- 
haviour is governed by the OPE of the 
currents through which the hard momenta 
flows. Again, this step is well defined in- 
dependently of the large-_/V c expansion; in 
practice, however, the calculations simplify 
when restricted to the appropriate order in 
the l/iV c -expansion one is interested in. In 
particular, chiral loops are subleading in the 
1 /-/V c -expansion. 

3. Introduce a large-N c approximation for the 
underlying Green's functions. 

As already mentioned, Green's functions 
in QCDqo are meromorphic functions. The 
minimal hadronic approximation (MHA) 
that we are proposing consists in limiting 
the number of poles to the minimum re- 
quired to satisfy the leading power fall-off 
at short-distances, as well as the appropri- 
ate xPT long-distance constraints. 

The three steps above can be done analytically, 
which helps to unravel the intricacies of the un- 
derlying dynamics. The method is, in princi- 
ple, improvable (unlike the other models dis- 
cussed above). It can be improved by adding 
more constraints from the next-to-leading short- 
distance inverse power behaviour in the OPE 
and/or higher orders in the chiral expansion. This 
has been tested within a toy model of QCD^ in 
ref. j§ 

We have checked this approach with the calcu- 
lation of a few low-energy observables: 

i) The electroweak Am, mass difference which 
we have already discussed JjJ . 

ii) The hadronic vacuum polarization contribu- 
tion to the anomalous magnetic moment of the 
muon a^. The MHA in this case requires one 
vector-state pole and a perturbative QCD contin- 
uum. The absence of dimension two operators in 
QCD in the chiral limit, constrains the threshold 
of the continuum. The result thus obtained |H| is 
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compatible with the more precise phenomenologi- 
cal determinations which use experimental input. 

Hi) The 7T° — » e + e~ and rj — > decay rates. 

These processes are governed by a (PVV) three- 
point function, with the Q 2 -momentum flowing 
through the two ^-currents. The MHA in this 
case requires a vector-pole and a double vector- 
pole. The predictions of the branching ratios 

^-frpz^y are in good agreement with the present 
experimental determinations. 

These successful predictions have encouraged 
some of us to pursue a systematic analysis of 
if-physics observables within the same large-iV c 
framework. So far, the following calculations have 
been made: 

• The Bk -Factor in the Chiral Limit p4[ . 

• Weak Matrix Elements of the Electroweak 
Penguvn Operators Q 7 and Q, H& 

• Hadronic Light-by-Light Scattering Contri- 
bution to the Muon g — 2 |36[|3^ ] . 

• Electroweak Hadronic Contributions to the 
Muon g- 2 |§. 

We also have preliminary results on: 

• The Long-Distance Contribution to the 
Kl — * fJ- + Decay Rate 39 . 

• Weak Matrix Elements of the Strong Pen- 
guin Operators Q4 and Q 6 EO]. 



5. FACTOR B K OF K° K° MIXING 

The factor in question is conventionally defined 
by the matrix element of the four-quark operator 
Qas=2(x) = {sLj tJ, d L )(sLl f j,dL)(x): 

(K a \QAS=2(0)\K°) = ^f 2 K M 2 K B K (p) . (20) 

To lowest order in the chiral expansion the oper- 
ator Qas=2(^) bosonizes into an 0(p 2 ) term: 

- ^5as =2 (m) [(D^)U] 23 [{D^)U] 23 , (21) 

4 V ' 

d» /f°a K o - ( r M) 23 ( r ) 23 + . . . 

with <7as=2(a0 the low energy constant (which 
depends on the renormalization scale ji) to be 



determined. A convenient choice of the underly- 
ing Green's function here is the four-point func- 
tion W^^ R (q,l) of two left-currents {shl^dh) 
and (slT^l), which carry the g-momentum one 
has to integrate over, and two right-currents 
{d~R"f a SR) and (c£r7^ ' sr) with soft /-momentum 
insertions, which couple to the external (r a )32— 
and (r^)32-sources. The coupling constant 
3AS=2(A t )> which has to be evaluated in the same 
renormalization scheme as the Wilson coefficient 
Cas=2(m) has been evaluated, is then given by 
the integral Q 



5as= 2 (m) = 1 



1 



32tt 2 F 2 



(4tt M 2 



t/2 



r(2- e /2)7 (Q 2 y/ 2 



dQ 2 



dn g \img^W£*'f R (q,l) 



W LL (Q 2 ) 



(22) 



MS 



conceptually similar to the one which determines 
the electroweak constant C in Eq. (fjj). The 
renormalization-scale invariant Bk~ factor is then 
given by the product 



Bk = ^C A s=2(m) x 5as=2(m) 



(23) 



The large— N c hadronic approximation of the 
Green's function Wll(Q 2 ) in Eq. (|2^), which 
fulfills the leading OPE short-distance con- 
straint and the long-distance constraints which 
fix Wll(0) and W LL (0) in x?T, requires one 
vector-pole, a double vector-pole and a triple 
vector-pole Q The integral in Eq. ( ^2|) can then 
be evaluated, with the result IM 



B K = 0.38 ±0.11 



(24) 



The /x-scale and renormalization scheme de- 
pendence in 3as=2(m) cancels with the one in 
Cas=2(a*)j when evaluated at the same approxi- 
mation in the //7V c -expansion. 

When comparing this result to other deter- 
minations, specially in Lattice QCD, it should 

4 In fact, this goes beyond the strict MHA which here only 
requires a vector-pole. It is the extra input of Wll (0) and 
W' L r (0), as known from xPT, which allows us to improve 
on the strict MHA. 
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be realized that the unfactorized contribution in 
Eq. ( p2| ) is the one in the chiral limit. It is pos- 
sible, in principle, to calculate chiral corrections 
within the same large- N c approach, but this has 
not yet been done. 

The result in Eq. (p4|) is compatible with the 
old current algebra prediction jll| which, to low- 
est order in ;\PT, relates the .Br— factor to the 
K + — > 7r + 7r decay rate. In fact, our calculation 
can be viewed as a successful prediction of the 
K + — > 7r + 7r decay rate ! 

As discussed in ref. |42[ |, the bosonization of 
the four-quark operator Qas=2 and the bosoniza- 
tion of the operator Q2 — Qi which generates 
AI = 1/2 transitions, are related to each other 
in the combined chiral limit and next -to-leading 
order in the l/7V c -expansion, when mixing with 
the penguin operators is neglected. Lowering the 
value of Bk from the factorized 0(A 2 ) predic- 
tion: Bk = 3/4, to the result in Eq. (p4|), is cor- 
related with an increase of the coupling constant 
<78 in the lowest order effective chiral Lagrangian 
(see Eq. (H)) which generates AI = 1/2 transi- 
tions, and provides a first step towards a quan- 
titative understanding of underlying dynamics of 
the AI = 1/2 rule. 

6. ELECTROWEAK PENGUINS 

We shall be next concerned with the four-quark 
operators generated by the so called electroweak 
penguin diagrams of the Electroweak Theory 



• C 7 {p)Q 7 + C 8 (^)Q S . 



(25) 



with Cj(fi), Cs(/j.) the Wilson coefficients of the 
operators 

Ql = Z{s L 'fd L ) e «few), (26) 

q— u,d,s 

Q 8 = -12 e q (s L q R )(q R d L ) . (27) 

q— u,d,s 

These operators generate terms of 0(p°) in the ef- 
fective chiral Lagrangian j4j|; therefore, their ma- 
trix elements, although suppressed by an e 2 fac- 
tor, are chirally enhanced. Furthermore, the Wil- 
son coefficient Cs has a large imaginary part in- 
duced by the top-quark integration, which makes 



the matrix elements of the Qg operator to be par- 
ticularly important in the evaluation of e'/e. 

Within the large-iV c framework, the bosoniza- 
tion of these operators produce matrix elements 
with the following counting 



(Qr)\o {p o) = Q(N c ) + 0(N°). 



(Q S )\o( P o) = 0(N 2 C 



O(N c 



(28) 



(29) 



Zweig suppressed 



6.1. Bosonization of Qj 

The bosonization of the Qj operator to 0(p°) 
in the chiral expansion and to 0(N C ) is very sim- 
ilar to the calculation of the ^-contribution to 
the coupling constant C in Eq. (Q). An eval- 
uation which also takes into account the renor- 
malization scheme dependence has been recently 
made in ref. |55| with the result: 

(Q 7 )b(pO) - 6 (0\(s L ^d L )(d R ^ SR )\0) x 

V v ' 

<Oi(m)> 

tr^AfV^y , and A^ 23) = 8 i2 S j3 . (30) 
Here, the vev (Oi(fi)) is given by the integral 
in ( \\ ( 3(6-3) (4V) e / 2 

(0l(/A) H 16, 2 r(2- £ /2) x 

r°° 

/ dQ 2 (Q 2 )^/ 2 (-q 2 u lr (q 2 ))} (31) 

with Hlr(Q 2 ) the same correlation function as in 
Eq. (||). In QCD^ this integral can be, formally, 
evaluated exactly. When restricted to the same 
MHA which has been used to calculate Am w , one 
gets the simple result 



(Oi 



32tt 2 



My 



,(32) 



where A 2 =/i 2 exp(l/3 + k), with k= —1/2 in the 
naive dimensional rcnormalization scheme (NDR) 
and k = +3/2 in the 't Hooft-Veltman scheme 
(HV). 
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6.2. Bosonization of Qg 

To lowest O(p ) in the chiral expansion, the 
four-quark operator Qg bosonizes as follows 

(Q8>|o( p o)=-12<0 2 (/i)>tr (u\^ 3) U^Q R )\(33) 

{Q\(s L s R )(d R d L )\0) 

As noted in refs. the vev 02^) also ap- 

pears in the Wilson coefficient of the 1/Q 6 term in 
the OPE of the same Hlr{Q 2 ) correlation func- 
tion as in Eq. (||), for which the MHA to QCD^ 
gives a rather good approximation, as we have 
already discussed. This offers the possibility of 
obtaining an estimate of the vev O2(a0 beyond 
the large- N c approximation, where C?2=> \(i'' l P) 2 , 
and, therefore, without having to fix a value 
for the (V'V') _corl densate, which is poorly known. 
Unfortunately, as one can judge from the results 
in Table |l|, different methods to extract the value 
of O2(a0 lead, at present, to rather different val- 
ues for the weak matrix elements: 

M 7 , 8 ee ((7T7r) z=2 |Q 7 , 8 |0) , at /i = 2 GcV. (34) 

In view of the diversity of results obtained with 
the dispersive approach, depending on which type 
of sum rules are used; and the unknown system- 
atic errors of the quenched approximation and the 
chiral limit extrapolations in the lattice results, 
it looks like it will take some time before one can 
claim that these matrix elements are known reli- 
ably. 

6.3. Test of the Zweig Rule 

Besides the important issue of getting an 
accurate determination of the Qs matrix el- 
ements, a reliable determination of the vev 
(0|(sLSfl)(dfl.G?z,)|0) would be most welcome as 
a test of the Zweig rule in the scalar and pseu- 
doscalar sectors. More precisely, 

(0 2 ( M )> = + (0\(s L s R )(d R d L )\0) c , (35) 

where the unfactorized contribution (the second 
term) involves Feynman diagrams which require 
gluon exchanges between at least two quark loops. 
These are the so called Zweig-suppressed contri- 
butions, which are indeed 0(N®) in the 1/N C 
expansion. There are reasons to suspect that 



subleading terms in the 1/N C expansion involv- 
ing Green's functions of scalar (and pseudoscalar) 
density operators might be important, unlike 
those which only involve vector (and axial- vector) 
currents. There are several phenomenological ex- 
amples of this: the rj' mass, the possible existence 
of a broad a meson, large final state interactions 
in states with J = and 1 = 0, etc. We are 
considering the possibility that the appropriate 
expansion for these exceptional Green's functions 
could be a 1/N C expansion in which rif/N c is 
held fixed, where nj denotes the number of light 
flavours; the kind of expansion originally advo- 
cated by G. Veneziano fl45| . 

Formally, the connected part of the vev 
{^\{sLSR){d R di J )\Q) c is defined by the integral 

with ^5>ij(Q 2 ) the two-point function 

l - J d 4 xe lq - x {<0|T [dd(x) 88(0)] |0) - 

(0|T [d 7B d(aO s l5 s(0)] |0)} . (37) 

It would be very helpful to get some information 
on the vev (0\(sLS R )(d R dL)\0} c from lattice QCD. 
Ultimately, this could provide a way to focus on 
the origin of the discrepancies in Table § 

7. GLUONIC PENGUINS 

Finally, I shall make a few analytic remarks 
concerning the sector of the four-quark opera- 
tors generated by the strong Penguins of the Elec- 
troweak Theory 

■■■C i ( t i)Q i + C 6 ( t i)Q 6 , (38) 

with Ci(p), Cq(^j) the Wilson coefficients of the 
operators 

04 = 4 ^ {s L ^q L ){q Ll ^d L ), (39) 

q— u,d,s 

Qe = S {s L q R ){q R d L ) ■ (40) 

q—u,d : s 

These operators generate contributions to the 
coupling constant gs of the 0(p 2 ) chiral La- 
grangian in Eq. (|^). The Wilson coefficient C*6(/i) 



10 



Table 1 

Matrix Elements Results for Mj^ (see Eq. (j34|)) in GeV 3 using different methods. 
METHOD M 7 (NDR) M 7 (HV) 



Ms (NDR) M 8 (HV) 



Large— iV c mha 
Knecht, Peris, de Rafael 



0.11 ±0.03 



0.67 ± 0.20 2.3 ±0.7 



(with aj corrects.) |35|]46 



2.5 ±0.. 



Dispersive Approach 

Narison^ 0.17 ±0.05 

Cirigliano et al g| 0.16 ± 0.10 

Bijnens et al @ 0.24 ± 0.03 

Cirigliano et al (OPAL data) |o) 0.21 ± 0.05 



0.49 ±0.07 
0.37 ±0.08 



1.4 ±0.3 
2.2 ±0.7 
1.2 ±0.8 
1.7 ±0.3 



2.5 ±0.7 
1.3 ±0.8 



Lattice QCD 
Bhattacharia et al ]5l[ | 
Donini et al ]52| 
RBC coll. || 
CP-PACS coll H 



(0.32 ±0.06) 
0.11 ±0.04 
(0.27 ±0.03) 
(0.24 ±0.03) 



0.18 ±0.06 



(1.2 ±0.2) 
0.51 ±0.10 
(1.1 ±0.2) 
(1.0 ±0.2) 



0.62 ±0.12 



Only the results obtained with methods which can exhibit an explicit dependence on the renormalization 
scale are quoted in this Table. The numbers in brackets have been obtained after informal private 
discussions with some of the authors of the various lattice collaboration. 



has also an important imaginary part generated 
by the integration of the heavy flavour degrees 
of freedom, which makes the operator Qq partic- 
ularly relevant for the evaluation of e'/e in the 
Standard Model. 

Seen from the large-A^ c framework point of 
view, the only terms which have been retained 
(so far) in analytic evaluations of weak matrix el- 
ements of four-quark operators, is best illustrated 
by showing the terms which, equivalently, would 
be retained in the lowest order anomalous dimen- 
sion matrix of these operators. Non-zero entries 
in this matrix appear then in three blocks, in the 
sector of the Qi, Q2, Qi, Q& and Qs operators, 
as follows: 





3 



3 

N, 





1 "/ 

3 N c 
1 "f 
3 N c 



1 "/ 
3 N c 

' 3 N c 



\ 



— iL_ , 



(41) 



More precisely, the mixing between matrix el- 
ements of the (Q 2 , Qi)-operators on the one 
hand and those of the penguin-like operators 
on the other, is neglected. Furthermore, in the 
{Q2, Qi)— sector, only the leading and next -to- 
leading terms in the l/iV c -expansion have been 
retained 



2 34 



Mixing between the strong pen- 
guin sector and the electroweak penguin sector 
is also neglected, but terms O(^f) in the mix- 
ing of Q4 and Qq are retained. Exceptionally, for 
reasons already discussed above, the Zweig sup- 
pressed O^/iV^-corrections in Q% are also re- 
tained. I can now discuss some recent analytic 
observations which have been made. 



7.1. Final State Interactions and e'/e 

The importance of final state interactions in 
the evaluation of K — > mr matrix elements has 
been reexamined in a series of recent papers |5j| ■ 
The main observation is that a simple numerical 
estimate of e'/e can be obtained if one proceeds 
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as follows: 

i) The evolution of the relevant Wilson coeffi- 
cients from the M^-scale to the m c -scale is re- 
tained, exactly, to next to leading order in pQCD. 

ii) The bosonization of the relevant Qe and Q% 
operators is done at the leading 0(N 2 ) only. No- 
tice that this corresponds to the further restric- 
tion where the matrix in Eq. |ll]) becomes diag- 
onal, with only non-zero entries in the (Qe, Qs) 

sector: ( ^ ^ V Starting at this level, the 

only 1/N C contributions which are then retained, 
are those generated by an approximate estimate 
of the the tttt final state interactions. 

iii) The tttt final state interactions are approx- 
imated by an Omnes-like resummation of the 
leading tttt chiral loops with non-zero disconti- 
nuities. The authors of ref. [|55[ justify this pro- 
cedure by the fact that it reproduces rather well 
some of the phenomenologically known 0(p 4 ) lo- 
cal terms of the electroweak chiral Lagrangian. 

The prediction thus obtained 



xlO" 3 , (42) 



where their estimates of the various sources of er- 
rors are indicated in underbracing, agrees within 
errors with the present experimental world aver- 
age |§ from NA31, NA48 and KTeV: 




Re(e//e) 



Exp. 



(1.66 ±0.16) x 10" 



(43) 



The obvious question here is to know how sta- 
ble is this prediction when terms of 0(N C ) are 
also retained in the bosonization of Q 6 and Q$. 
The scenario with large deviations from the naive 
0(N 2 ) bosonization of Qq and Qg, which cancel 
in their overall contribution to e'/e, cannot be 
excluded. Let us not forget that the restriction 
to 0(N 2 ) terms in the bosonization of the four- 
quark operators, fails dramatically to reproduce 
the AI = 1/2 rule; a fact which, perhaps, should 
not be so surprising since, after all, it is only when 
the 0(N C ) terms are incorporated as well, that 
the QCDoo-properties of planarity apply. 



7.2. Bosonization of Qq 

Related to the previous discussion is the ques- 
tion of the bosonization of the Q6~operator. It 
has been known for sometime that this opera- 
tor, to 0(N 2 ), gives a contribution to the cou- 
pling constant g$ in Eq. (^) which is modu- 
lated by the product of the ratio (iJip) 2 / 'F§ times 
the 0(p 4 ) L5 -coupling of the strong chiral La- 
grangian (known from the f-nl Ik ratio). It can be 
shown that, in the presence of unfactorized terms, 
this contribution is modified as follows fl4Gl : 



9b\q 6 =C&{p) 



A) 



47T/i 2 



e/2 



2tt 2 F 4 L(2 - e/2) 



dQ 2 (Q 



2\l-e/2 



Wdg(Q 2 ) 



MS 



(44) 



where the function Wdg(Q 2 ) is the analog of 
the function W L l(Q 2 ) in Eq. (||) and the 
function Il LR (Q 2 ) in Eqs. (@) and (|). Here, 
Wdg{Q 2 ) is the invariant function associated 
with a four-point function VV^g^^q, /) of two 
density-currents {sl^r) and (qnclr,), which carry 
the ^-momentum one has to integrate over, and 
two right-currents (d]ij a un) and (ur^^sr) with 
soft Z-momentum insertions, which couple to the 
external (r Q )i 2 - and (rg ) 3:L -sources. We have 
found that the effect of the unfactorized piece in 
Eq. ( |44| ) is very important. The large effect can 
already be seen in the low Q 2 behaviour of the 



function Wdg{Q 2 ) in X?T (B 



Jim W DG (Q 2 ) 



l f 



(45) 



Q 2 

where the factor L 5 — |L 3 ~ 11 x 10~ 3 , which gov- 
erns the behaviour of Q 2 Yv 'dg(Q 2 ) at the origin 
(the integral of the first term vanishes in dimen- 
sional regularization) , is very large as compared 
to the individual values of L5 and L3. This large 
value results mostly from the effect of the low- 
est vector state which the first 0{N 2 ) term in 
Eq. ([l5]) does not see. Clearly, this has serious 
implications both for the AI =1/2 rule and e'/e. 
An evaluation, along these lines, of the coupling 
g s \ Qs is in progress. 
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8. CONCLUSIONS and OUTLOOK 

We think that large-./V c QCD provides a very 
useful framework to formulate approximate cal- 
culations of the low-energy constants of the ef- 
fective chiral Lagrangian, both in the strong and 
electroweak sector. 

The approach that we propose has been tested 
successfully with the calculation of a few low- 
energy observables discussed in section [|. Here, 
we have described various applications to the 
evaluation of weak matrix elements which, for the 
first time in analytic calculations of this type, ex- 
plicitly show the cancellation of the renormaliza- 
tion scale between the short-distance contribu- 
tions and the long-distance contributions. The 
analytic implementation of the approach is suf- 
ficiently simple so as to provide an explanation 
of the dominant underlying physical effects. It 
may also help, as a guidance, to unravel impor- 
tant physical effects in Lattice QCD numerical 
simulations. 
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